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Abstract 

First, a condition for the Bose-Einstein transition in the superfluid 
hquid hehum is presented due to the creation of a free neutron spin- 
less pairs in a liquid helium and a dilute neutron gas mixture. We 
proposed a new model of dilute Bose gas, where presence atoms in 
the condensate is a suppressor for the collective modes as well as a 
creator for single-particle excitations On other hand, it is shown that 
the terms, of the interaction between the excitations of the Bose gas 
and the density modes of the neutron, meditate an attractive interac- 
tion via the neutron modes, which in turn leads to a bound state on a 
spinless neutron pair. The lambda transition is defined by a condition 
for the Bose-Einstein transition, which transforms refiected neutron 
pair modes to single neutron modes. 
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1. Introduction. 



The motivation for our theoretical study of the very low-temperature 
properties of the dilute hard sphere Bose gas is an attempt at a microscopic 
understanding of superfluidity in helium ^He. We proceed by discussing some 
experimental and theoretical investigations. 

The connection between the ideal Bose gas and superfluidity in helium 
was first made by London [1] in 1938. The ideal Bose gas undergoes a 
phase transition at sufficiently low temperatures to a condition in which the 
zero-momentum quantum state is occupied by a finite fraction of the atoms. 
This momentum-condensed phase was postulated by London to represent the 
superfluid component of liquid ^He. With this hypothesis, the beginnings of 
a two- fluid hydrodynamic model of superfluids was developed by Landau [2] 
where he predicted the notation of a collective excitations so- called phonons 
and rotons. 

The purely microscopic theory with mostly utilizes-technique, was first 
described by Bogoliubov [3] within the model of weakly non-ideal Bosc-gas, 
with the inter-particle S- wave scattering. Based on the application of the 
presence of a macroscopic number of condensate atoms Nq N (where N is 
the total number of atoms), Bogoliubov has dropped a density operator term 
which describes the fluctuation of atoms above the zero momentum level, 
the Bogoliubov obtained the dispersion curve for single particle Bogoliubov 
excitations (Bogoliubov phonon-roton modes). 

The dispersion curve of an excitations excited in superfluid helium has 
been accurately measured a function from momentum [4]. Within this exper- 
iment, in lambda transition, the position of a sharp peak inelastic neutron 
scattering intensity defines by energy of the single particle excitations, and 
there is appearing a broad component in inelastic neutron scattering inten- 
sity, at higher momenta of atoms. For explanation of the appearance of a 
broad component in inelastic neutron scattering intensity, the authors of pa- 
pers [5-7] proposed to consider the presence of collective modes in the super- 
fiuid liquid ^He, which are represented a density excitations. In this respect, 
the authors of this letter predicted that the collective modes represent as the 
density quasiparticles [8]. In these works, these density excitations and den- 
sity quasiparticles are appeared due to remained density operator term for 
describing atoms above the condensate, which was neglected by Bogoliubov 
[3]. 

The last time, the authors of letter [9] discovered an existence of scattering 
between atoms of the superfluid liquid helium, at lambda transition, which is 
confirmed by calculation of the dependence of the critical temperature on the 
interaction parameter as scattering length. On other hand, as it is indicated 
by authors of [4], there are two type excitations in superfiuid helium, at 
lambda transition. These facts imply that it needs to revise the concept 
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of determination condition for Bose-Einstein condensation in the superfluid 
liquid helium. Obviously, the peak inelastic neutron scattering intensity is 
connected with a registration of neutron modes by neutron-spectrometer, 
which in turn defines sort of excitations. Therefore, in this letter, we present 
to determine a condition for Bose condensation, at lambda transition, by 
a registration neutron-spectrometer a single neutron modes or neutron pair 
modes. 

In this letter, we present a new model of a nonideal Bosc gas for describing 
of the superfluid liquid helium. The given model is based on the application 
of the Penrose-Onsager definition of the Bosc condensation [10] which is 
based on the condition for a condensed fraction of atoms ^ = const. The 
later explains a broken of Bosc-symmctry law for the atoms of the Bosc 
gas in the condensate level, and gives a fully exact solution to the model 
of dilute non-ideal Bose gas, at quantity of the condensate fraction < 
^ < 1. In this context, the new model of a nonideal Bose gas, presented 
herein, leads to an absence of the collective modes because appearance of 
atoms in the condensate is the suppressor for collective modes as well as the 
creator for single-particle excitations. Therefore, we prove that the density 
excitations and the density quasiparticles, proposed by authors of letters 
[5-8], are unphysical. On other hand, at selection condition for condensed 
fraction ^ ^ 1, we may state that new model of a dilute Bose gas describes 
a thermodynamic property of a realistic superfluid liquid helium. 

Further, we investigate a helium liquid-dilute neutron gas mixture where 
exists the term of interaction between the Bogoliubov modes and the density 
modes of the neutrons, which due to application a canonical transformation 
of the Hamiltonian of system, the term of the interaction between the density 
of the Bogoliubov modes and the density of the neutron modes is removed 
by meditated an effective attractive interaction between the neutron modes, 
which in turn determines a bound state on neutron pair. 

2. New model of a dilute Bose gas. 

For beginning, we present a new model of a dilute Bose gas for describing 
property of superfluid liquid helium. The given model considers a system 
of N identical interacting atoms via S-wave scattering. These atoms, as 
spinless Bose-particles, have a mass m which are confined in a box of volume 
V . The main part of the Hamiltonian of such system is expressed in the 
second quantization form as: 

1 

= E «P + ^ E UpQpQ'^ ( 1) 

Here cit and fip- are, respectively, the "creation" and "annihilation" op- 
erators of a free atoms with momentum p; is the Fourier transform of a 
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S-wave pseudopotential in the momentum space: 



U,= '-^ (2) 

m 

where d is the scattering amphtude; and the Fourier component of the 
density operator presents as 

Pi 

According to the Bogohubov's theory [3], it is a necessary to separate the 
atoms in the condensate from those atoms fiUing states above the condensate. 
In this respect, the operators ciq and are replaced by c-numbers ciq = 
^0 — within approximation the presence of a macroscopic number of 
condensate atoms A^o ^ 1- This assumption leads to a broken of Bose- 
symmetry law for atoms occupying in the condensate state. To be refused a 
broken from Bosc-symmctry law for bosons in the condensate, we apply the 
Penrose-Onsager's definition of the Bose condensation [10]: 

No 

lim — = const (4) 

No,N-^oo N 

This reasoning is a very important factor for microscopic investigation of 
the model non-ideal Bose gas because the presence of a macroscopic number 
of atoms in the condensate means Nq » Np^o (where Ap-^o is the occupation 
number of atoms with momentum p^O). In this respect, we may postulate a 
following approximation for an occupation number of atoms with momentum 
p: 

No-*oo iVg 

The next step is to find the property of operators , by applying 
(5). Obviously, 

ai -. 

JV0-.00 ^ ' 

and 

limj^^ = <5,„, (7) 

Excluding the term pi = 0, :the density operators of bosons and gi 
take the following forms: 

gp^^o{a1:p+V2cp) (8) 
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and 



= 0Vo(a-P + V2ct) (9) 



where and ct are, respectively, the Bose-operators of density-quasiparticles 
presented in reference [8] which in turn are the Bose-operators of bosons used 
in expressions (6) and (7): : 

'^P~ fnw S ^pi-p^Pi ~ ~~fK X/ '^PiiP^Pi ~ ~~fK^P (1^) 



and 



Thus, we reach to the density operators of atoms and ^t, presented 
by Bogohubov [3], at approximation ^ = const: 

Qp= ^o(atp+a^^ (12) 

and 

g±^^o(a-p+a±) (13) 



which displays that the density quasiparticles are absent. 

The identical picture is observed in the case of the density excitations 
which was predicted by the Glyde, Griffin and Stirling [5-7] where was pro- 
posed a presentation gp in a following form: 



Qp = \JNq [aZp +ap+ Qpj ( 14) 

where terms involving pi 7^ and , Pi 7^ p are written separately; and the 
operator describes the density-excitations: 

" ~7w ^ (15) 

At inserting of (6) and (7) into (15), the term, representing as the density- 
excitations, vanishes because = 0. 

Consequently, the Hamiltonian of system, presented in (1) by support of 
(12) and (13), reproduces an extension form of the Bogoliubov Hamiltonian, 
at approximation ^ = const: 

Ha^Y.[7T- + 'nfiv^] O'fap + ^ E ( ^-p^f + ^p^-p) (1^) 
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where v = y = \J '^'^'^l-^y^ is the velocity of sound in the Bose gas 
which depends on the density atoms in the condensate 

For evolution of the energy level it is a necessary to diagonalize the Hamil- 
tonian Ha which is accomplished by introduction of the Bose-operators &t 
and 6p> by using of the Bogoliubov linear transformation [3] : 

where is the unknown real symmetrical function of a momentum p. 
Substitution of (17) into (16) leads to 

p 

hence we infer that bt and hp are the " creation" and " annihilation" op- 
erators of a Bogoliubov quasiparticles with energy: 



, , ,^ . (19) 
\2mJ J ^ ' 

In this context, the real symmetrical function Lp of a momentum p is 



found 



It is well known, the strong interaction between the helium atoms is very 
important and reduces the condensate fraction to 10 percent or ^ = 0.1 
[4], at absolute zero. However, as we suggest a new model of dilute Bose 
gas, proposed herein, may have an significant application for describing of 
thermodynamic properties of the superfiuid liquid helium because the S-wave 
scattering between two atoms, with coordinates ri and r2 in the space of 
coordinate, is presented by the repulsive potential delta-function [/^ = ^""'^^ 
from r = ri — f2. On other hand, the presented model works on the condensed 
fraction ^ <C 1 in differ from the Bogoliubov model where ^ ^1. We can 
see that in statistical equilibrium, the equation for the density of charged 
bosons in the condensate is presented in a manner analogous to the one 
obtained by Bogoliubov [3]: 

^-^-lV-^-lvi±5W- (21) 



where b^bp is the average number of the Bogoliubov quasiparticles with the 
momentum p: 
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ekT — 1 

3. The effective attractive potential interaction between neu- 
tron modes in a helium liquid-dilute neutron gas mixture. 

We now attempt to describe the thermodynamic property of a hehum 
hquid-dilute neutron gas mixture. In this context, we consider a neutron 
gas as an ideal Fermi gas consisting of n free neutrons with mass m„ which 
interact with A'^ interacting atoms of a superfluid hquid helium. The helium- 
neutron mixture is confined in a box of volume V. The Hamiltonian of a 
considering system Ha^n consists of the term of the Hamiltonian of Bogoliubov 
excitations Ha in (18) and the term of the Hamiltonian of a dilute Fermi 
neutron gas as well as the term of interaction between the density of the 
Bogoliubov excitations and the density of the neutron modes: 

Ha,n + ^J^fip +1^77^2 UQQpQ-p,n (22) 



with 



P 2^ P^O 



p,(7 " p 



where ai^ and a^^a are, respectively, the operators of creation and annihi- 
lation for free neutron with momentum p, by the value of its spin z-component 
a —ll |; Uo and Ui are, respectively, the Fourier transforms of the repulsive 
interaction between the density of the Bogoliubov excitations and the density 
modes of the neutrons as well as between neutron modes: 



and 
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where do is the scattering amplitude between a helium atoms and neu- 
trons; di is the scattering amplitude between neutron modes; /j, = is 
the relative mass. 

Hence, we note that the Fermi operators dt^ and d^^^ satisfy to the Fermi 
commutation relations [•••]+ as: 



+ P,P 



■ Ka' (24) 
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(25) 



[^p>,a+, ,]+ = (26) 
The density operator of neutrons with spin a in momentum p is defined 

as 

Qp,n = J2 ^i-p,a^Pi,<y (27) 

where gt^ = g.p^n 

The operator of total number of neutrons is 

p,a 

On other hand, the density operator, in the term of the Bogohubov quasi- 
particles Qp included in (22), is expressed by following form, to application 
(17) into (12): 

Hence, we note that the Bose- operator hp commutates with the Fermi 
operator ap^„ because the Bogoliubov excitations and neutrons are an inde- 
pendent. 

Now, inserting of a value of operator ^^-from (28) into (22), which in turn 
leads to reducing the Hamiltonian of system Ha,n'- 



Ha,n — X/ 2rYl ^^'^'^P''^ ~^ 2V ^ ^^Qp,nQ-p,n + 
p,a p 

+ Y.^P^k+^^Y.J^X^^-p + k)o^^^^^ (29) 

p p 1 p 

Hence, we note that the Hamiltonian of system Ha,n iii (29) is a similar 
to the Hamiltonian of system an electron gas-phonon gas mixture which was 
proposed by the Prolich at solving of the problem superconductivity (please, 
see the Equation (16) in H. Frolich, Proc.Roy Soc, A215, 291-291 (1952) 
in the reference [10] ), contains a subtle error in the term of the interaction 
between the density of phonon modes and the density of electron modes 
which represents a third term in right side of Equation (16) in [10] because 
the later is described by two sums, one from which goes by the wave vector 
w but other sum goes by the wave vector k. This fact contradicts to the 
definition of the density operator of the electron modes g^; (please, see the 
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Equation (12) in [11]) which in turn already contains the sum by the wave 
vector k, and therefore, it is not a necessary to take into account so-called 
twice summations from k and w for describing of the term of the interaction 
between the density of phonon modes and the density of electron modes Thus, 
in the case of the Prolich, the sum must be taken only by wave vector w, due 
to definition of the density operator of electron modes with the momentum 
of phonon w. 

To allocate anomalous term in the Hamiltonian of system Ha,n, which 
denotes by third term in right side in (29), we apply the Prolich approach 
[10] which allows to do a canonical transformation for the operator Ha^n 
within introducing a new operator H: 

H = exp (^+) Ha,n exp (30) 
which is decayed by following terms: 



H = exp (^S'+) Ha,n exp = Ha,n - [S, Ha,n] + \ [^, [^, Kn]] " " " " (31) 

where the operators represent as: 

^+ = E4 (32) 
p 

and 

S = T.Sp (33) 
p 

and satisfy to a condition S'^ — —S 
In this respect, we assume that 

4— ^p(^pA--^^n^p-) (34) 

where Af is the unknown real symmetrical function from a momentum p. 
In this context, at application Sf from (34) to (33) with taking into account 
Q-p,n = Qp,n: ^^en we obtain 

= E 4 = E ^P^P> (b=P - bf) (35) 
p p 

In analogy manner, at gt^^n — 6p,n^ we have 

= E 4 = E ^P^Sn - K) = - E ^pQp,n (b-P - b}) (36) 
p p p 
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Obviously, the given form of the operator S in (35) coincides with a 
presentation of one, which was included by the Frdlich in Equation (18) by 
reference [10], within the function ip{k, w) from momenta k and w, which in 
turn represents as the c-number and equals the Lp{k, w) = v?(w), where ip{w) 
is the real function from w which is a similar to the value proposed in this 
work. The later result is connected with a presentation of a correct form, for 
the term of interaction between the density phonon modes and the density of 
electrons, connected with the Frolich Hamiltonian in [10], as it is mentioned 
in above. 

To find A^, we substitute (29), (35) and (36) into (31). Then, 



p V p p 

1, 



j[S, [S, Ha,n]] = E A'iepg^,ng-p,n (38) 
p 

and [S, [S, [S, -^a,n]]] = within application a Bose commutation relations 

S-S \Qpi,m Qp2,n\ ~ and [cij^_^^^Clj^-^^^u , Qp2,n\ ~ 0. 

Thus, the form of new operator H in (31) takes a following form: 



p p V » V / 



- T^Yl ^pU i/iVo W Y^T^ ^P." ^-P>" ~ E ^P^p{i-p + Q-p,n + 
p V p p 

+ J2 Alepg^^rig-p,n + E ^^h^P''' + ^ UiQp,nQ-p,n (39) 
p p,cr p 

The transformation of the term of the interaction between the density of 
the Bogoliubov modes and the density neutron modes is made by removing 
of a second and fifth terms in right side of (39) which leads to obtaining of a 
quantity for Ap: 



A.-'-^-Jl^ (40, 



2epV V 1 - ^P 



In this respect, we reach to reducing of the new Hamiltonian of system 
(39): 



H = E^p*?V+E^«^-«p> + ^E^i^p>^-p>- 

p p,<T p 
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- ^ H ^P^O Y^^Jp,ne-p,n + ^i^pQp,nQ-p,n (41) 

p VP p 

As result, the new form of Hamiltonian system takes a following fom: 

H = '£s^±hp+Hn (42) 
p 

where Hn is the effective Hamiftonian of a neutron gas which contains an 
effective interaction between neutron modes: 

Hn^Y. ^^S<^^P> + ^ E ( + ^1 ) Qp,nQ-p,n (43) 

p,a " p 

where is the effective potential of the interaction between neutron 
modes which takes a following form at substituting a value of Ap from (40) 
into (41): 



Vp = -2ApU^^J \^ + 2A^^pV = '-^ ^ (44) 

V 1 - Vep(l - Lpj 

In this letter, we consider following cases: 1. At low momenta atoms of 

a helium p << 2mv, the Bogoliunovs quasiparticles in (19) represent as the 

1— ^ 

phonons with energy ep ~ pv which in turn defines a value Lp i+T. ~ 

^1 — in (20) or Lp ^ 1 — In this context, the effective potential 

between neutron modes takes a following form: 



2 mU^No _ Annuel 
Vp^ p 

The value ei is the effective charge, at a small momenta of atoms: 



yp---^---^ (45) 



2. At high momenta atoms of a helium p » 2'mv, we obtain ep ~ 
^ + mv'^ in (19) which in turn defines L^- pa in (20). Then, the effective 
potential between neutron modes presents as: 



Vp'^ p'^ 
where 62 is the effective charge, at high momenta of atoms: 



^2 = 7^ 



2h\ Vt: 
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Consequently, in both following form presents the effective scat- 

tering between two neutrons in momentum space: 

= ^ (47) 

where e* = ei, at small momenta of atoms; and e* = 62, at high momenta. 

4. Creation Spinless Neutron Pairs and Single Neutrons. 

We now find the potential interaction between two neutrons in the coor- 
dinate space: 

= hj:(yp + Ui)-e''^ = ^ E(^- + ^i)-^^''' (48) 

where = is defined by wave number w. 

Hence, we take into consideration a condition for S-wave scattering be- 
tween neutron modes within using of a condition Wfdi <^ 1 (where Wf — 

^3^^ is the Fermi wave number for neutron gas; di is the scattering ampli- 
tude between neutrons). This reasoning implies that the scattering between 
two neutrons is presented in the coordinate space by a following form: 

Vif^ - ^ E • e-" = 47r r (V^ + U,)w''-^^^dw (49) 
K ~r Jo \ J wr 

where we introduce a following approximation as ^''^^'"^ ~ 1—^^^^ because 

the conditions w < Wf and Wfd <^ 1 lead to Wfr <^ 1 {wf — ^^y^^ is 
the Fermi wave number). The later condition defines a state for distance r 



V 



between two neighboring electrons f <^ 

In this context, by taking into consideration = we obtain 

2 2/'37r^n\3 Uin f Trel n Ui [ n\o ^ ^ 



V(r> . -Ifa^e; ^— J +^ + - ^ ^-J J .r' (50) 

This approximation means that there is an appearance of a screening 
character in the effective scattering because as we see the later depends on 
the density electron modes. 

We now attempt to investigation of the effective Hamiltonian of a neutron 
gas in (43) is rewrite down in the space of coordinate by following form: 

- 2 

Hn^±H,^ --^ E A. + E ni n - f, I) (51) 

i=l i=l i<3 
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where Hi is the Hamiltonian of system consisting two neutrons with op- 
posite spin which have a coordinates and fj: 



The transformation of considering coordinate system to the relative coordi- 
nate f—fi—fj and the coordinate of center mass R — ^-^y^i ^® have 

Hi = -^^R - —K + V{r) (53) 

To find the binding energy E < oi neutron pair, we search the solution 
of Schrodinger equation with introduction of wave function ip{r): 

H.ijsif) = Eijj.ir) 

In this respect, we have a following equation 

Ari^sir) + V{f)il;s{r) = ^^(r) (54) 

which may determine the binding energy £' < of neutron pair Inserting 
value of V{r) from (50), and denoting E — Eg 



Ar - 

. m„. 



/37r^n\3 Uifi 



*\ V 



2V 



El 

2 



We transform the form (55) by a following form: 



+ A-6'VU,(r) = 



where we take 



/nir. 



7re? n 



h \ 3 V 



and 



Hi 
2 

V 



Um 



(55) 



(56) 



By application of the wave function '0s (r) via the Chebishev-Hermit func- 
tion Hs{it) from an imaginary number as argument it [12] (where i is the 
imaginary one; t is the real number; s = 0; 1; 2; • • •), the equation (61) has a 
following solution as: 
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where 



at ^ < within 



dt' 



Consequently, the quantity of the binding energy Es of neutron pair with 
mass mo = 2mn is rewritten as: 



/rrir, 



\ ?, V 2\V 



1 



2„2 



s + - -167r^e: 



at s = 0;l;2;••• 
The normal state of neutron pair corresponds to quantity s = which 
defines the maximal quantity of the binding energy of neutron pair: 



Eq 



h 



7re? n 



Ui(n\3 2 2 



V 



Um 
W 



< 



(58) 



which implies that the creation of neutron pair is appeared by the following 
condition: 



o o /Stt n\ 3 Uin , , 

> ^ (59) 

Thus, the spinless neutron pair with mass mo = 2m„ is created in a 
helium liquid-dilute neutron gas mixture, at condition (59), which determines 
a bound state on a neutron pair with binding energy (58), The later depends 
on the density of atoms in the condensate and therefore, may define the 
temperature of lambda transition T^, at condition 



3 2 \v) 
At low momenta atoms of a helium p << 2mv, 



(60) 



V 327r2[/o2 m VF. 
but at high momenta atoms of a helium p » 2mv, 



(61) 



TVo^T, ^ 3f/i r ( n 



V l&T:W^m\V 



(62) 
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Thus, in the state of temperatures < T < Tx, there are a creation 
neutron pairs but at T > T^, the neutron pair is decayed on two free neu- 
trons because the condition (59) is broken because the condition (59) takes 
a following form 



which leads to positive meaning of Eq > 0. 

In conclusion, we note that the new model of Bosc gas, presented in 
this letter, might be useful for describing of the thermodynamic properties 
of a dilute gas of the Boson-Fermion mixtures confined in traps. Hence, 
we note that the correction form of the Prolich Hamiltonian system, for 
describing property of a phonon gas-electron gas mixture in [13], may lead to 
creation spinlcss electron pairs, in addition of a phonon gas. The presented 
model of a helium liquid-dilute neutron gas mixture gives explanation of the 
presence of a broad component in inelastic neutron scattering intensity, at 
lambda transition [4] because there is an appearance single neutron modes, in 
addition to neutron pair modes, which in turn determine a superfluid phase 
for helium. 




(63) 
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